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ABSTRACT 

Numerical modeling of the Pulsed Inductive Thruster exercising the magnetohydrodynamics code, 
MACH2 aims to provide bilateral validation of the thruster’s measured performance and the code’s 
capability of capturing the pertinent physical processes. Computed impulse values for helium and argon 
propellants demonstrate excellent correlation to the experimental data for a range of energy levels and 
propellant-mass values. The effects of the vacuum tank wall and mass-injection scheme were investigated 
to show trivial changes in the overall performance. An idealized model for these energy levels and 
propellants deduces that the energy expended to the internal energy modes and plasma dissipation 
processes is independent of the propellant type, mass and energy level. 


INTRODUCTION 

The Pulsed Inductive Thruster (PIT) is a unique propulsion system due to a combination of its distinctive 
physical operation, demonstrated performance and potential advantages over all other electric rockets. 1 
A spiral induction coil is powered by a series of capacitors in Marx-loop configurations and produces a 
strong azimuthal electric field once the pulsed current is passed through. The field breaks down the 
injected propellant over the surface of the fiat coil and the associated radial magnetic field induces an 
azimuthal current within a thin gas layer (figure 1). This current interacts with the magnetic field to 
produce an axially-directed electromagnetic force. Such inductive acceleration circumvents the need for 
electrodes and the erosion-related lifetime limitations that plague traditional electromagnetic thrusters. 
The PIT can operate with any propellant, however single-shot operation demonstrated optimum 
performance with polyatomic -molecule fuel. Specifically, the PIT-MkV operating with ammonia at 4.6kJ 
demonstrated nearly-constant efficiencies exceeding 50% for a wide range of specific-impulse values; 
4000s<Isp<8000s. Its propensity to work best with such fuels renders the thruster ideal for water- 
propellant operation and associated missions such as In-Situ Propellant Utilization (ISRU) where 
refueling can be realized by appropriate water-rich cometary rendezvous or by capitalizing on Europa’s 
oceans for deeper space assignments. Its pulsed operation can provide elevated maneuverability — the 
aforementioned operation produced an impulse range of 0.05N-s<I-bit<0.12 N-s - which in conjunction to 
its water-propellant affinity establishes the concept very attractive for Orbital Express Architectures 
during which satellites can refuel and upgrade electronics by docking to on-orbit modules. 

In order to provide added confidence in the PIT’s demonstrated capabilities the magnetohydrodynamics 
(MHD) code, MACH2 is employed to model its operation and offer preliminary validation by 
comparisons with experiments. Specifically, the paper aims to validate the code’s capabilities in capturing 
the pertinent physics by comparing experimental results to the numerical findings. 
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Figure 1. — The Pulsed Inductive Thruster (PIT) (Courtesy TRW) and 
a schematic of the relative acceleration process due to the 
dominant interaction of the applied radial field 
and the induced azimuthal current. 


NUMERICAL MODEL 

MACH2 is a time-dependent, two-dimensional, axisymmetric, multi-material code that can be applied to 
problems of complex geometries due to its multi-block structure. 2 The computational mesh can move in 
an Arbitrary-Lagrangian-Eulerian (ALE) fashion allowing applicability to both diffusive-and 
dispersive-dominated problems as well as code validation. The mesh can be refined via a variety of 
adaptive schemes to capture regions of varying characteristic scale. The set of the single-fluid, MHD 
equations is time-advanced with finite-volume spatial differencing, and the boundary conditions are 
applied via the ghost-cell technique so that no special conditional statement is necessary at the 
boundaries. 

The mass continuity and momentum equations assume a compressible, viscous fluid with the latter 
including both real and artificial viscosity effects. The stress-deviator can be chosen to evolve under 
elastic stress for strength of material calculations 3 or modeled as a Newtonian fluid to upgrade the code to 
a Navier-Stokes solver. The code includes two ablation models that allow mass addition due to solid 
evaporation and have been successfully employed to model ablation-fed pulsed plasma thrusters. 4 The 
electrons, ions and radiation field are in thermal non-equilibrium, so MACH2 solves three energy 
equations. These include thermal conduction with anisotropic transport 5 and three different models for 
radiation cooling. 6 Evolution of the magnetic field is prescribed by the induction equation that includes 
resistive diffusion, the Hall effect and the thermal source for magnetic fields. Various models for the 
plasma resistivity are available. They comprise classical anisotropic resistivity, 5 several anomalous 
resistivity models and contributions from electron-neutral collisions applicable to weakly ionized gases. 2 
In many engineering applications the source of magnetic flux is applied currents produced from 
externally-applied voltage differentials. For this, the code includes a variety of circuit models such as 
LRC, Pulse-Forming-Networks, sine-waveforms and several others. Other additional physical models 
include laser-pulse energy deposition and a detonation package. The set of the MHD equations is 
completed by an equation of state that can be either analytic or tabular. The latter is provided by the 
SESAME library that includes semi-empirical models for the thermodynamic properties, transport 
coefficients, (including opacities) and average ionization state under local thermodynamic equilibrium. 
These models have been constructed and are being maintained by the T-l and T-4 groups at Los Alamos 
National Laboratory. 7 
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The level of sophistication and capability of the MACH2 code has been instrumental in providing 
invaluable insights to a variety of plasma problems. Some of these include plasma opening switches, 8 
inertial-confinement fusion and alternative concepts, 9 compact toroid formation and acceleration 
schemes, 10 gas and solid density z-pinch implosion physics, 11 laser-target interactions, 12 high-power 
plasma source design, 13 magnetic nozzles 14 and a variety of plasma thrusters. 15,16 Its diverse success 
establishes the code as a primary numerical tool toward the understanding of the intricate physical 
processes and feasibility issues of the proposed propulsion concept. 


PIT PHYSICAL MODELING 

A substantial body of experimental data is available for the Pulsed Inductive Thruster (PIT) comprising a 
diverse range of propellants, energy levels and propellant mass values. 1 The power supply consists of a 
series of 18 capacitors in Marx-loop pairs that are charged in parallel. This results in an effective 
discharge voltage that is double that of the each capacitor’s charge voltage which in turn provides the 
necessary high electric fields for breakdown. In particular, the configuration results in an effective 
capacitance and voltage that relate to the single capacitance and charge voltage through the total energy 
available: 


E o =^( #of capacitors )CF 0 2 =^C eff v} ff 


( 1 ) 


where V 0 =l/2V e ff. For the polyatomic propellant experiments the effective capacitance was 9(liF while the 
earlier efforts with monatomic propellants utilized smaller capacitors with effective capacitance of 4.5qF. 
The resulting oscillatory current waveform provided pulse duration of the order of 15psec with peak 
currents exceeding lOOkA for a range of charge voltage of 10kV<F o <16kV and propellant mass of 
1.5mg<m<10mg. Mass is introduced through a nozzle (figure 2) by a fast opening gas valve. For optimal 
distribution of the injected gas the mass-pulse is short enough such that it ceases as the leading edge 
reaches the confining lexan. The PIT’s performance data exhibited approximately-constant-efficiency 
trends for a wide range of specific impulse values (500s<Isp<9000s) providing impulses exceeding 
0.1 N-s. Monatomic propellants demonstrated efficiencies on the order of 20% while polyatomic 
propellants with the higher-capacitance bank dramatically improved exhibiting maximum efficiencies 
above 50% with ammonia. 

The MACH2 simulations utilized the helium and argon data for comparisons as the SESAME tables 
provided a complete real equation of state that includes a degree of ionization semi-empirical model. The 
physical model included thermal non-equilibrium of a single-viscous fluid with classical transport. 
Boundary conditions modeled thermal and magnetic field insulators with no-slip. The computational grid 
(figure 2) extended well downstream of the thruster’s exhaust region to comprehensively capture the 
acceleration process and assure no influence of the outlet boundary conditions that model variables at 
zero gradient. Grid resolution was maximized in the vicinity of the coil to assure capture of the fast-rising 
field’s diffusion and associated gradient. In particular, the significant acceleration process occurs during 
the rise-time, t r ~0.9qsec, which implies a characteristic diffusion depth of about 1.78cm for a 2eV 
plasma, (S~[rp r ] /2 , where rj-1000/7^ 3/2 m 2 /s is the electrical diffusivity with T e in eV). The axial grid-cell 
dimension, Az=1.56mm assures accurate capture of the field’s gradient with more than 10 cells resolving 
the diffusion depth. 
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Figure 2. — Schematic of the half-plane thruster with the 
computational grid utilized, (to scale). 1.) Nozzle with 
pulsed mass-valve, 2.) Conical pylon, 

3.) Confining cuff (lexan). 


Table 1. — Comparison of relevant characteristics of 
the PIT and LRC emulated waveforms. 


Vef fi E: 

J maxi trise • 

First zero, t Q : 

J min • 


PIT 

28kV, 1764J 
lOOkA, 0.75ps 
3ps 
-40kA 


MACH2-LRC 
28kV, 1764J 
lOOkA, 0.9|lis 
2.7|lis 
-22.6kA 


Initial conditions assigned uniform density in the vicinity of the coil (the confined region shown by the 
highest-resolution grid in figure 2) to at best resemble the evolved gas injected from the valve-and 
uniform room temperature. The significance of non-uniformities present due to the injection scheme is 
addressed later. This implies that ionization of the gas which occurs after the breakdown is calculated by 
the MACH2 ionization model consistently with the evolution of the rest of the pertinent variables. 
Emulation of the current waveform used the LRC circuit model. Specifically, a circuit model had been 
developed 1 to calculate the PIT’s current waveform with significant accuracy. 15 For the MACH2 
simulation the LRC external parameters were adjusted to as-best match the important features of the 
aforementioned model. The comparison is shown in table 1 and the LRC current waveform used to 
simulate helium operation is depicted in figure 3. 

It should be emphasized that the PIT circuit model included a semi-empirical statement of the plasma’s 
inductance axial variation along with the plasma’s resistance influence. The MACH2 LRC model even 
though it does include the influence of the latter, it does not presently include the inductance variation. 

Thrust calculation was performed at the outlet boundaries as the integral of the dynamic pressure and 
integrated to provide the impulse. This implied that even though the effective pulse duration was about 
8ps computational times needed to be in excess of lOOps so as to capture exhaust of all the propellant 
through the downstream boundaries. The thrust histogram for the representative case of V () = 1 4kV and 
m=3mg is shown in figure 4 along with the evolution of the impulse. The sharp rise of the thrust value at 
about 30psec suggests exhaust of a compact plasma consistent with a slug-acceleration notion. 
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Figure 3. — Calculated Current waveform used 
to emulate the PIT’s circuitry, F^28kV. 


Figure 4. — Calculated thrust histogram and 
associated impulse evolution, V e ^=28kV, m=3mg. 


This was confirmed by inspection of the two-dimensional density distribution as calculated by the code; 
they depicted a fairly sharp mass outline that included the majority of the propellant. The second peak in 
the evolution proposes additional acceleration of part of the remaining mass by the second peak at 
minimum current. The reduced gradient especially after the peak suggests a more dispersed plasma that 
takes longer to pass through the exhausting boundary. The width of the thrust pulse also indicates the 
degree of effective coupling and of course the eventual exhaust speed. In other words, a narrower 
histogram with higher peak thrust values would indicate better coupling and thus increased exhaust 
velocities for the same propellant mass and energy level. In turn, the aforementioned observations imply a 
substantially effective coupling of plasma and magnetic flux that allows induction of significant currents 
and consequently electromagnetic force. Such behavior is quite encouraging as it suggests that MACH2 
captures the relevant coupling mechanisms as expected from such concept. 


COMPARISONS WITH EXPERIMENT 

Computations using the MACH2 code were performed for the range of available charge voltage and 
propellant mass. Specifically, experiments addressed a range of energy levels between 900J and 1764J 
with propellant mass variation from 0.75mg to 9.2mg. For each pair, data for the impulse were provided 
which in turn produced the range of specific impulse and efficiency values. MACH2 was invoked to 
compute all the different arrangements for a full range comparison of Impulse versus specific energy 
illustrated in figures 5 through 8. The computed impulse correlates very well with the experimental 
values. MACH2 computed impulse values for the 900J case exhibit an overestimation of about 10% when 
compared to the single experimental point. This discrepancy is partially attributed to an experimentally 
observed critical-mass value below which incomplete breakdown was observed along with the consequent 
reduction in impulse. Even though the code can not address breakdown physics, this implicit incomplete 
coupling of plasma and magnetic flux is partially reflected by the computations as indeed efficiency is 
decreasing with decreasing propellant mass value. However, the magnitude of this reduction is not fully 
captured by the code’s physical models and thus the overestimation. It is of interest to note though, that 
MACH2 ? s partial capture of this critical-mass trend is indicative of additional influencing processes other 
than breakdown behavior. 
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Propellant mass (mg) 

Figure 5. — MACH2 impulse comparison with 
experimental data for helium and 1296J. 



Propellant mass (mg) 

Figure 6. — MACF12 impulse comparisons with 
experimental data for argon and 1764J. 



Propellant mass (mg) 

Figure 7. — MACF12 impulse comparison with 
experimental data for helium and 900J. 
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Figure 8. — MACH2 impulse comparisons with 
experimental data for helium and 1764J. 


Even though the comparisons with experimental data are extremely encouraging certain simplifications 
and deficiencies of the code at its present state may be offsetting so as to result in the computed accuracy. 
In particular, three issues arise. Experimental data were obtained in a restricted vacuum tank of 1.2m in 
diameter (the thruster’s diameter is lm) while the simulations implemented a free (outlet) boundary 
condition for the exhausting plasma. (See figure 2, the right outermost boundary of the computational 
region.) The simulations assumed uniform initial conditions over the surface of the induction coil-which 
may tend to overestimate the calculated impulse values as uniform flux coupling is encouraged-and the 
lack of the LRC circuit model reaction to the plasma dynamics resulted in a current waveform that 
underestimates the second current peak by almost half. The latter two issues may be actually offsetting 
each other as any underestimation of the acceleration from the second current peak may be equalized by 
the aforementioned inaccuracy due to uniform initial conditions. For this the first two issues were 
addressed. 


A new set of simulations was contrived that would emulate the vacuum tank’s walls at 0.1m away from 
the confining lexan. This is accomplished by changing the right boundary of the computational region to a 
wall, no slip boundary condition. The effects and setup are clearly shown by the two-dimensional 
distributions of the velocity vectors (figure 9) where the latter depicts the simulation with the confining 
vacuum-tank wall. A comparison of the computed performance is shown in table 2 and implies 
insignificant influence from the restrictive vacuum tank. 
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Figure 9. — Velocity vector distribution depicting the 
influence of the wall boundary condition. 

Top: Outlet, Bottom: Emulation of 
vacuum-tank wall at the 
right boundary. 

Table 2. — Calculated impulse values with and 
without influence from the vacuum-tank wall. 


EfJVmfmg) 

900/4.2 

1296/3 

1296/1.5 

1764/3 


I (outlet) 
44.3 mN-s 
35.7 mN-s 
26.5 mN-s 
42.2 mN-s 


I (tank wall) 

43.6 mN-s 

35.7 mN-s 

28.7 mN-s 
43.1 mN-s 


The second issue that can be addressed by the code in its present state is the mass injection scheme. The 
thruster utilizes a gas valve and a conical nozzle to guide the cold propellant gas toward the induction 
coil. The duration of the pulse in the experiments was designed such that the desired amount is injected 
before the leading edge of the expanding gas reaches the confining cuff. A delay before discharge 
followed such that all the propellant is given the opportunity to distribute itself over the induction coil. 
Emulation of the process by the MACH2 code introduces the mass at the conical nozzle as helium at 
room temperature and subsonic speeds. This choice of the mass-flow-rate value ensured injection of the 
total propellant in 2msec. (The particular simulation considered the 1.5mg case.) The two-dimensional 
mass density distributions depicted in figures 10 and 11 outline the evolution and behavior of the cold 
propellant as it expands toward the induction coil. At 2msec the flow rate is interrupted and the emulated 
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discharge is initiated which is somewhat dissimilar from the experimental setup that enforced a delay 
such that the mass is distributed as uniformly as possible. In other words, the MACH2 simulation will 
address non-uniformities from the mass -injection scheme that are much more profound than during actual 
operation. One of the main differences is that the MACH2 scheme prevents a significant portion of the 
mass to interact with the fast-rising magnetic field, (the propellant still flowing along the pylon in 
figure 11, 2msec). 

Comparisons of the calculated impulse values between the uniformly distributed properties and this 
injection scheme for V () = 1 2k V and m=1.5mg show insignificant influence of the latter. Specifically, 
uniform initial conditions predict an impulse of 26.5mN-s while the more pragmatic mass -injection 
scheme calculates 23.9mN-s. The discrepancy is insignificant, especially considering that during the 
experiments the mass was given the opportunity to distribute itself as uniformly as possible before the 
voltage drop was applied. A plausible explanation arises from interrogation of the axial distributions of 
the magnetic field temperature and ionization level, but can as well be deduced by noting that only a 
portion of the propellant (of the order of the diffusion depth) adjacent to the coil participates in the 
magnetic flux coupling and current induction. This hotter, ionized fraction is subject to electromagnetic 
acceleration, but as it propels downstream it entrains the remainder colder propellant. This colder 
propellant’s location is relatively inconsequential (whether it resides immediately downstream of the 
accelerating fraction or still expanding over the conical pylon) as it will be eventually swept by the fast 
moving portion. 

The aforementioned deductions can be extrapolated as to the effect of the LRC model’s plasma- 
inductance deficiency which resulted in underestimating the minimum current’s magnitude by half. They 
imply that any additional acceleration due to the second current peak is minimal. If indeed the 
overestimation from assuming uniform initial conditions was originally equalized by the deficiency in 
calculating the acceleration due to second current peak then the latter’s contribution is only of the order of 
(26.5-23.9)mN-s=2.6mN-s~10%. This is minimal of course partly due to the lower energy levels involved 
in the helium experiments and thus the relative significance of the plasma’s inductance. The LRC model 
in MACH2 needs to be amended to include the effects of plasma dynamics before higher-energy-level 
modeling is performed. 



Figure 10. — Mass-density evolution as it is Figure 11. — Mass-density evolution as it is injected 

injected from the conical nozzle, from the conical nozzle, (1.2msec and 2msec). At 

(0.3msec and 0.8msec). 2msec the emulated discharge is initiated. 


NASA/CR — 2003-2 12714 


8 



THEORETICAL MODELING 


In order to gain further insights into the acceleration mechanism of the Pulsed Inductive Thruster (PIT), 
we will construct an idealized model by assuming a purely radial magnetic field that is radially and 
azimuthally uniform. This can be formulated by the following expression (see figure 12), and boundary 
conditions: 


B r {z)= ^° J f{l) where C, = — and /(o) = l;/(l) = 0;^ 7 /(C)|^_i = /'(l) = 0 (2) 

r 0 ~ r i z o 

The function, y(Q is arbitrary provided it satisfies the stated boundary conditions subject to an effective 
current conduction zone, z 0 . Application of Ampere’s law produces an expression for the induced 
azimuthal current, 


- 1 dB 

VxB = - r 


-VoJ'e =>h 


jf%) 


z „ 3? """ h-n'ko 

Then the purely axial electromagnetic force exerted on the plasma can be formulated: 

1 2 " '• n„J 2 AQft) 


= j j»B r dV = J J j- 


V 


0 0 r 7 - 


Vo~ r i) Vo~ r i ) o 


(3) 


(4) 



Figure 12. — Schematic of the Pulsed Inductive 
Thruster (PIT) utilized for the idealized model. 


The value of the integral involving the axial distribution is independent of the function, fiQ so long it 
adheres to the boundary conditions; it is -1/2. Hence, the accelerating force is, 

F{f)= ^0 4£fc_t5) = I rj( , )2 = 

- 'i I 2 2 (r 0 -r,) 

where V is the inductance per unit length. 
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If we further assume a LRC circuit configuration for the PIT we can write the governing differential 
equation as follows: 


V(t) - L e j(t ) - R e J(t) - V p (0 = 0; (?) = L p (?) j{t) + [R p (?) + L p (?)]/(?) (6) 


where the subscript, e denotes the circuit’s external elements and subscript,/? denotes plasma parameters. 
For the idealized approach we will assume constant effective values for the plasma elements. We can 
identify the plasma's impedance as Z = R p + L p , and since the plasma’s inductance is L p = L'z 0 we can 
combine to produce the idealized circuit equation: 


V(t)-Lj(t)-[R e +Z\j(t) = 0 where L = L e +L'z 0 
The solution to the differential equation is the familiar LRC circuit waveform, 


(7) 


*7(/) = "~ ex p| 
to L 


\-{R e + z )~ 

sin(co?); (0 = ,- 

2 L 

V 


LC 


ALr 


( 8 ) 


that allows calculation of the following convenient integral to be utilized throughout the remainder of the 
model’s formulation, 


|j(?) 2 J? = j(^exp 


(Re+Z) . 

2 L 


sin(co^) 


CV 2 

dt = 2 — - 


J 


0 0 

Then the total impulse generated by the accelerator is given by 


2{R e +Z) (R e +Z) 


(9) 


I = \F Z (? )dt = U'j J{tfdt = , (10) 

0 0 \ e ) 

and since I=mU e we can formulate an expression for the plasma’s kinetic energy, KE and the energy 
dissipated to the circuit’s external elements via resistive loss, E R 

KE= X -mUl= {L )2e2 ° ; E R =R e ]j{t) 2 dt = j (11) 
2 8 m(R e +Z) 2 J 0 ( R e + Z ) 

A simple statement of conservation of energy that distinguishes these two energy modes can be written as 
E R +KE=(l-QE 0 , where ^ is the fraction of the total energy deposited into the plasma via resistive heating, 
and dissipated due to irreversible processes. Substitution of the expressions stated above results in a 
second-order polynomial with respect to the quantity, Re+Z, the positive root of which can be easily 
expressed as follows: 
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( 12 ) 


(1 -^\R e +Z) 2 -R e {R e +Z)- 


L' 2 E ( 

8m 


— = 0 => R„ + Z = ■ 


R„ + a/ R. 


il 


+ {\-^)L' 2 E 0 /2m 


2 ( 1-0 


Substitution into the impulse expression with some algebraic manipulation results 

/ = (l -QL'E 0 Jbn . L , = K\X ( Xr 0 +r i ) 

R e y[2m + ^2mi?g +(1 + 0 ^ /2 £' o 2 ^° ~ ^ ^ 


(13) 


where the expression for the inductance per unit length is repeated for convenience. It is important to 
emphasize that the above manipulations imply the fraction ^ to be independent of the plasma’s 
impedance. We also note that as the external circuit element losses are diminished (R e ->0) or as the 
propellant mass, m is decreased, impulse tends toward a maximum that is a function of the square root of 
the propellant mass and is independent of the plasma’s inductance per unit length. The impulse expression 
allows calculation of the specific impulse, Isp=I/mg and efficiency=I2/2mEo along with verification by 
comparisons to experimental data. The latter is accomplished by utilizing the experimentally deduced 
external resistancel, 7? e =5mQ, along with the PIT’s geometry (r o =0.5m, ri=0.2m) for a range of propellant 
masses and energy levels operating with different propellants. Such comparisons are displayed in figures 
1 3 through 1 5 for helium, argon and carbon dioxide propellants which show impulse versus initial energy 
normalized over propellant mass. 



Specific Energy, E 0 /m (J/mg) 


Figure 13. — Comparison of the idealized model 


with experimental data for helium. ^=0.78. 



Figure 15. — Comparisons of the idealized model 
with experimental data for carbon dioxide. ^=0.78. 



Figure 14. — Comparisons of the idealized model 
with experimental data for argon. ^=0.78. 



Figure 16. — Efficiency versus specific impulse 
for argon propellant at £’ 0 =1764J. 
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The comparisons for this energy range highlight a quite unexpected trend. The fraction of the total 
available energy expended in heating the propellant through resistive heating (i.e., deposition of energy to 
heat and ionize the gas) and all other dissipating phenomena (i.e., heat transfer to the surroundings, 
viscous momentum losses, etc.) is independent of the propellant mass, energy level or even propellant 
type, it is constant at 0.78! In other words, the idealized model insinuates that for all these propellants the 
efficiency can not exceed 22% regardless of propellant mass value and energy level. In addition, the 
impulse expression indicates that the efficiency is not constant, but rather approaches a constant 
maximum value as propellant mass is decreased. This is due to the fact that the effect of external losses 
(due to R e ) diminishes as propellant mass is decreased. It is instructional to plot efficiency versus specific 
impulse which is displayed in figure 16 for argon propellant. We note the model’s trend as it 
appropriately predicts that in the absence of any available kinetic energy, efficiency vanishes. Further 
development of the model is necessary in order to decipher the behavior of this energy dissipation factor 
at higher energy levels, but most importantly for ammonia and hydrazine that have produced much higher 
efficiency values than 0.22. 


CONCLUSIONS 

The magnetohydrodynamics code, MACH2 was employed to model the performance of the Pulsed 
Inductive Thruster (PIT) operated with helium and argon propellants within a range of 900J-1764J and 
propellant mass range of 0.75mg to 9.2mg. Comparisons between predicted and experimentally measured 
impulse values demonstrated excellent agreement. The effects of a restrictive vacuum tank used to obtain 
the experimental data were examined by simulating a wall, no-slip boundary condition at the location of 
the vacuum tank. Performance was insignificantly affected. Effects of the mass injection-scheme were 
also investigated for one energy and mass case to also show that performance was minimally affected by 
the non-uniform distribution of the injected propellant over the induction coil. An idealized model 
produced an expression for the impulse based on a fraction of energy expended to the internal plasma 
energy modes and dissipation effects. Comparisons with experimental data for this energy range 
(i? o <2000J) imply that this energy fraction is invariable regardless of propellant type (for three different 
propellants; He, Ar and C0 2 ), propellant mass and energy level. 
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